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ABSTRACT 



We report on an investigation of various problems related to the theory 
of the electroweak phase transition. This includes a determination of the 
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nature of the phase transition, a discussion of the possible role of higher 
order radiative corrections and the theory of the formation and evolution of 
the bubbles of the new phase. We find in particular that no dangerous linear 
terms appear in the effective potential. However, the strength of the first 
order phase transition is 2/3 times less than what follows from the one- loop 
approximation. This rules out baryogenesis in the minimal version of the 
electroweak theory. 
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1. With the recognition that baryon number violation is unsuppressed 
at high temperature in the standard model has come the realization that the 
electroweak phase transition might be the origin of the observed asymmetry 
between matter and antimatter p], ^ |3| . In order to have sufficient departure 
from equilibrium, it is necessary that this transition be rather strongly first 
order. As a result, there has been renewed interest in understanding un- 
der what circumstances the transition is first order, and how the transition 
proceeds. The minimal standard model almost certainly cannot produce the 
observed asymmetry: it has too little CP violation, and, as we will see, its 
phase transition is too weakly first order for a Higgs more massive than the 
present experimentalphalimit. Nevertheless, for considering the features of 
the phase transition, it is a useful prototype, because it is weakly coupled 
and comparatively simple. 

Despite its apparent simplicity, understanding the phase transition in this 
theory has turned out to be surprisingly complicated, and the literature now 
contains contradictory claims on almost every point. In the first papers on 
this subject it was assumed that the phase transition is second order [|]. 
Later it was shown that if the Higgs mass is sufficiently small, the phase 
transition becomes first order ||. The problem of bubble formation was 
considered in some detail; it was argued that at least initially the bubble 
walls typically are rather thick @. While bubble wall evolution was only 
touched upon in these early efforts 0, simple arguments suggested that, e.g. 
for a 50 GeV Higgs, the motion of the wall would be non-relativistic. 

Recently, quite different views have been expressed about all of these 
issues. Brahm and Hsu, in an interesting paper [[7] , have argued that infrared 
effects spoil the one-loop analysis and claim to reliably establish that the 
transition is second order. Anderson and Hall || have thoroughly considered 
a number of aspects of the phase transition. Most of their results are in 
qualitative agreement with the earlier treatments 0, [J and with our previous 
results |J, while providing some quantitative improvement. However, they 
argue that the initial bubbles are thin. Gleiser and Kolb [ID"|, and Tetradis 



E| argue that fluctuations are so large that the transition does not proceed 
through the formation of (critical) bubbles. There has also been controversy 



as to whether or not bubble wall motion is ultrarelativistic f|l3l IHI . 
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In the present note, we will attempt to deal with these various issues. 
We will focus on the minimal standard model. As we have said, this model 
is not realistic, but we expect that the arguments and methods described 
here can be extended to more realistic situations. We focus on this case only 
because it is the simplest - and even here, we will frequently have to content 
ourselves with rather crude calculations. We will be able to give at least 
partial answers to each of the questions raised above. We will argue that, 
for sufficiently small coupling, the phase transition can be reliably shown to 
be first order. In particular, the linear term in the potential found in Ref. 
[0] is not present. On the other hand, we will see that the transition is more 
weakly first order than suggested by the one loop analysis. We will show that 
the bubbles are thick when they form, and the bubble nucleation rate cannot 
be reliably computed in the thin wall approximation. We will show that, for 
the range of Higgs masses considered here, the transition does in fact proceed 
by nucleation of critical bubbles. Finally, we will make some estimates of the 
bubble wall thickness and velocity. We will see that the early estimates of 
Ref. || are only reliable if particle mean free paths are very long. In practice, 
the relevant mean free paths are short, and the velocity of the bubble wall 
somewhat larger than these early estimates suggest. The problem of bubble 
wall propagation turns out to be surprisingly complicated, and we mention 
some of the issues which must be considered. We will see, however, that, 
for the first order transitions being considered here, the motion of the wall 
tends to be non-relativistic. In this paper, we will outline our treatment of 
each of these issues and describe the major results; details will be given in a 



subsequent publication JTJ]. 



2. We first consider the question of the order of the phase transition. 
The standard approach to this problem consists of computing the effective 
potential to one loop order. At high temperature, the one-loop expression 
for Vt is given, to a good approximation, by 

V{<f>,T) = D(T 2 - T 2 )0 2 - ET<f + . (1) 

Here A = m 2 H /2v 2 , and 

D = -^(2m 2 w + m| + 2m 2 ) , E = -^(2m^ + m|) ~ 1(T 2 , (2) 
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Here v = 246 GeV is the value of the scalar field at the minimum of V(<f>, 0) 
and lna B = 2\nAn - 2 7 ~ 3.51, \na F = 21nvr - 2 7 ~ 1.14 |L g §. 



The behavior of V(0, T) is reviewed in Refs. [ IB] , |T7| . At very high 
temperatures the only minimum of V(4>, T) is at = 0. A second minimum, 
0i, appears at T = Ti, where 

T2 - T o M - 3ET 

11 ~ 1 - 9Ey8X Tl D ' 01 " 2A T • (5j 

The values of V(<f>, T) in the two minima become equal to each other at the 
temperature T c , where 

Tc = i - E yx T D ' 0C = ^- (6) 

The minimum of V((f), T) at = disappears at the temperature T a , when 
the field in the second minimum is O = 3ET/Xt- 

3. The first order character of the transition is due to the 3 term 
in the potential. The appearance of such a term non-analytic in |0| 2 is the 
signal of an infrared problem, and raises concerns about the validity of the 
perturbation expansion. Indeed, these issues were raised in the early work 
on the subject [[|, [18], 19 1. The problem of infrared divergences at high tem- 



peratures in theories with light or massless particles has received extensive 
attention in the literature and is discussed in many textbooks fL6| , f20| . It 
is well known that these problems arise from the zero-frequency terms in 
the discrete frequency sums which appear in the computation of equilibrium 
quantities at finite temperature. The corresponding Feynman diagrams are 
those appropriate to a three dimensional field theory. For simplicity, we will 
consider here the contribution of W bosons. In the present case, the problem 
is most easily analyzed in Coulomb gauge, V • W — 0. The relevant propa- 
gators are those for the Coulomb lines, D 0Q , and transverse lines, D^. For 
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zero frequency, 



A» = i, D l3 (k)= 1 P^k) , (7) 
k k 2 + mwVP) 

where mw = #0/2, and = 5ij — The 'Goldstone' field and the 

physical Higgs scalar have standard scalar propagators with mass terms which 
are independent of the gauge coupling. The cubic term is readily extracted 
from the zero-frequency piece of the determinant. 2/3 of it arises from the 
transverse gauge bosons; the other 1/3 is obtained from the Coulomb line. 

In this gauge it is rather easy to see how higher orders of perturbation 
theory behave. At one loop, it is well-known that the Coulomb line, even 
for = 0, acquires a mass m 2 D = (N g + N c )g 2 T 2 /3 fl2~I| . In order to obtain 
a sensible perturbation theory for small 0, it is necessary to partially resum 
the perturbation expansion, i.e. to use in each order of the loop expansion 
a Coulomb propagator with this effective mass. In the infra-red, this resum- 
mation corresponds to integrating out all heavy modes (to one loop order), 
leaving only the u = modes, with a Debye mass. In practice, the tree level 
gauge boson mass at the minimum of the potential is small compared tomo. 
As a result, repeating the one loop calculation with this effective mass, the 
contribution to the cubic term from the Coulomb line disappears. The same 
is not true, however, for the transverse bosons. For zero 0, gauge invariance 
forbids a one-loop mass for these fields; the transverse polarization tensor is 

in fact given by Ily = 159 Pij- As a result, the contribution to the cubic 
term from these fields survives, and the net effect is to reduce the coefficient 
of the cubic term in eq. ([!]) by 2/3: Q 

E = -^—(2m 3 w + m 3 z ) . (8) 



Qttv 



This reduction means that, for a given Higgs mass, the phase transition 
is significantly less first order than one expects from the one loop analysis. 

1 After we obtained this result, we received a very interesting paper by Carrington ]2^| 
where the modification of the cubic term by high order corrections was also considered. 
Even though the author did not claim that these corrections reduce the cubic term by a 
factor of 2/3, after some algebra one can check that his result is equivalent to ours. 
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It has been pointed out that a minimal requirement of the phase transition 
is that the sphaleron rate after the transition be sufficiently small that the 
baryon number not be washed out. Using the (unimproved) one loop result, 
this gives a limit of about 42 — 55 GeV |23], PJ, at best just barely consistent 



with the present experimental constraints p4| . Allowing for the correction 
obtained here, the limit on m H is reduced by about 25%, clearly ruling out 
baryogenesis in the model. 

More generally, however, we can ask about the behavior of the pertur- 
bation expansion, particularly at small 0. Before making general remarks, 
it is helpful to consider the two loop diagrams involving transverse gauge 
bosons and scalars (Fig. 1). The zero frequency pieces of these diagrams, 
separately, give contributions ~ g s \<p\T 3 . If one combines these diagrams, 
however, being careful about combinatorics, they have the structure of an 
insertion of a polarization tensor on the transverse gauge boson line. Be- 
cause, as mentioned above, this tensor vanishes at zero momentum, the sum 
of these diagrams is less singular at small <f), and simply gives a correction to 
the quadratic term. We have checked all other diagrams at two loop order 
and shown that there are no linear terms in the potential. 

The authors of Ref. [[?]] found a linear contribution to the potential by 
simply substituting the scalar mass found at one loop back into the one loop 
calculation. Such a procedure is generally reliable when calculating Green's 
functions or tadpoles. Indeed, it is well known that the sum of the geometric 
progression, which appears after the insertion of an arbitrary number of 
polarization operators II(T) into the propagator (k 2 + m 2 ) -1 , simply gives 
(k 2 + m 2 + n(T)) _1 . However, this trick does not work for the closed loop 
diagram for the effective potential, which contains ln(k 2 + m 2 ). A naive 
substitution of the effective mass squared m 2 + II (T) instead of m 2 into 
ln(k 2 + m 2 ) corresponds to a wrong counting of higher order corrections. 
The simplest way to avoid the ambiguity is to calculate tadpole diagrams for 
^ instead of the vacuum loops, and then integrate the result with respect to 
cf). One can easily check by this method as well, that no linear terms appear 
in the expression for V(<p, T). 

The absence of linear terms does not automatically mean that higher 
order corrections are completely under control. A general investigation of 
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the infrared problem in the non-Abelian gauge theories at finite temperature 
suggests that the results which we obtained are reliable for 0^|T ~ T/3 



H| [19). Thus, a more detailed investigation is needed to study behavior 
of the theories with m#^10 2 GeV near the critical temperature, since the 
scalar field, which appears at the moment of the phase transition in these 
theories, is very small (see Fig. 2). However, we expect that our results are 
reliable for strongly first order phase transitions with 0^T, which is quite 
sufficient to study (or to rule out) baryogenesis in the electroweak theory. 

4. We turn now to the problem of bubble formation. At high tempera- 
tures, this becomes a problem in classical thermodynamics. One looks for a 
stationary point of the free energy, with the property that <fi —>■ as r —>■ oo, 
i.e. a solution of the classical field equations with potential Vq + Vr- Some 
care is required in solving this equation, however, since necessarily one is 
constructing a saddle point of the action (unstable modes corresponding to 
bubble growth or collapse). As a result, if one makes a poor approximation, 
one overestimates the probability of bubble formation. We believe this is the 
case of the analysis of Ref. ||, where formation of bubbles was studied in the 
thin wall approximation. This approximation works well if the height of the 
barrier between the two minima is much larger than the difference between 
the values of the effective potential in each of them. This is not the case for 
the phase transition with m# ^ 60 GeV. Numerical solution of the equations 
yields an action typically a few times larger than that obtained from the thin 
wall approximation. 

Even though the thin wall approximation fails, one can still study bubble 
formation analytically in a wide class of theories. Indeed, we have found 
that in the vicinity of the phase transition one can write, to a very good 
approximation, 

S 3 _ 38.8 D 3 ' 2 /AT\ 3 / 2 / 2 \ D AT \ 

T - E 2 ' \ T ) X v E 2 T /" ^ ' 

where 

r, n , «r 2.4 0.26 i . , 

/( a = 1 + - 1 + + . 10 

4 L I- a (i_ a ) z J 

Parameter a changes from to 1 in the range of temperatures for which the 
phase transition is possible. We have found that eq. ([TOD is correct with an 
accuracy about 1% in the most interesting range < a < 0.95. 
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5. In the case that the transition is weakly first order, it is natural to ask 
whether the transition actually proceeds through formation of bubbles, or if 
other sorts of fluctuations might be more important. In most treatments, it 
is assumed that the transition occurs once the bubble nucleation rate is large 
enough that the universe can fill with bubbles. In practice, because of the 
extremely slow expansion rate at the time of the transition, this means that 
the barrier is still high enough that the naive calculation of the nucleation 
rate gives an extremely small result; the three dimensional action is of order 
130 — 140. Given that the rate of formation of critical bubbles is so small, one 
might expect that other types of fluctuations which might equilibrate the two 
phases would be extremely rare. In Refs. |K| [LJ], [12], however, it has been 
argued that this is not the case. Roughly speaking, these authors arrive at 
this conclusion by estimating the mean square fluctuation of the scalar field 
about the symmetric minimum, 2 ms =< <f) 2 >, and comparing this with the 
value of the field at the other minimum. A rough estimate leads them to 
the conclusion that the < 2 >~ mT ~ 2 , so that it is not meaningful to 
consider the system as sitting in one vacuum or the other. Here m is the 
Higgs field mass near = 0. Subcritical bubbles, they argue, equilibrate the 
two phases even before one reaches the temperature T c . 

While we believe that for some range of parameters subcritical bubbles 
may be important, we do not believe that this is the case for the Higgs 
masses under consideration here. In estimating 4> rms , one should be careful 
to consider only long wavelength modes. Short wavelength modes will be 
associated with configurations with large gradient terms, which will collapse 
in a microscopic time. Also, one must be very careful with factors of tc and 2. 
A more detailed investigation based on the stochastic approach to tunneling 
gives the estimate for the amplitude of relevant fluctuations < 2 >~ 
[p5fl . Combining this estimate with our results for the mass m, i.e. for the 



curvature of the effective potential near = in the relevant temperature 
interval, T Q < T < T\, one obtains an estimate <fi rm ,s ~ -1 T. 

From this analysis, we see that the typical amplitude of the relevant scalar 
field fluctuations is substantially less than the separation of the two minima, 
unless the phase transition is very weakly first order. Even for rajj ~ 60 
GeV, the distance between the two minima remains five times greater than 
(firms- Including fluctuations with k 3> k max will give larger rms , but these 
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will collapse in a microscopic time, and will not serve to equilibrate the two 
phases. 



6. Understanding the motion of the wall is important mainly in the con- 
text of baryogenesis. We have already seen that no baryons will be generated 
in the single Higgs theory, unless the Higgs mass is well below the present 
experimental limit. We expect, however, that, for mu ~ 35 GeV (for which 
4>/T > 1), the bubble wall motion in this model will have many features in 
common with more realistic theories of baryogenesis, which require the phase 
transition to be strongly first order. Even in this simple model, determining 
the wall velocity and shape is a difficult problem, and we will content our- 
selves with rather crude estimates. In our analysis, we will assume that the 
wall achieves a steady state after some time. In particular, there is a frame, 
which we refer to as the 'wall frame', in which the scalar field and the parti- 
cle distributions are independent of time. We will assume that the principle 
source of damping of the walls motion is elastic scattering of particles from 
the wall; if there are additional sources of damping, they can only slow the 
wall even further. We will treat the velocity of the wall as a small parameter. 
We will have to check its validity a posteriori. With these assumptions, it is 
helpful to consider two limiting cases, depending on whether the size of the 
wall, 5, is large or small compared to the relevant mean free paths for elastic 
scattering, (. As an estimate of these mean free paths, we follow Ref. 
and take the longitudinal and transverse gluon propagators to include a mass 
proportional to m D . For top quarks, this yields a t = lQTta 2 s /3m 2 el . Multiply- 
ing by the flux one obtains ( ~ 4T _1 for quarks. For PF's and Z's the result 
is about three times larger. These numbers are consistent with results which 



have been obtained for the stopping power ||26|| . To get some feeling for the 
wall size, consider the system at temperature T c . At this temperature, the 
two phases can coexist, separated by a static domain wall. The <fi field in 
this domain wall is readily obtained by quadrature; for a 35 GeV Higgs, one 
finds 5-40 T-\ 

This suggests that the thin wall limit may not be a good one, but it is 
instructive to consider it in any case. In this limit, a typical particle passes 
through the wall, or is reflected from it, without scattering. An estimate in 
this limit was given in Ref. |J. If the problem is treated semiclassically, 
it is straightforward to calculate the extra, velocity-dependent force on the 
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wall, assuming that all particles approaching the wall from either side are 
described by an equilibrium distribution at some temperature. (Note that 
particles moving away from the wall are not described by an equilibrium 
distribution in this case; these particles are assumed to be equilibrated far 
from the wall, at some possibly different temperature and velocity). In this 
model it is straightforward to calculate the force on the wall to linear order 
in v. For bosons, the leading term is of order m 3 T, while for fermions, it is 
of order m 4 . Equating this force to the pressure difference on the two sides 
of the wall gives the velocity. Defining 

T c -T 

(11) 



T —T 



we can obtain an approximate equation for the velocity, valid for small e: 

•~£(TT3 <*"**>• <12) 

where the correction C comes predominantly from top quarks and is roughly 
of order one. For m t = 120 GeV and ran = 35 GeV (for which e ~ 1/4), we 
find v ~ 0.05. 

In this discussion we have made a variety of oversimplifications. One 
which is potentially important is our assumption of equilibrium densities on 
both sides of the wall. Some fraction of incoming particles, however, are 
reflected from the wall, and this will tend to lead to an enhancement of the 
particle density in front of the wall. For definiteness, consider top quarks. 
As a crude estimate, suppose the fraction of reflected particles is / (of order 
m/T), and suppose that the mean free path for processes which can change 
top quark number is r, with ( <C r. Then the density in front of the wall 
is enhanced by an amount of order fvn^J^. For the thin wall case, this is 
likely to reduce the velocity of the wall somewhat. In the thicker wall cases 
considered below, however, this is likely to be more important. 

Consider now the case 5 ^> (. In this case, to model the deviations from 
equilibrium due to the finite velocity of the wall, we can break the wall into 
segments of length (, and repeat the thin wall analysis for each of these 
segments. In particular, we assume that the distributions on either side of 
each segment are at equilibrium. We also assume that the temperature and 
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velocity of the particle distributions are constant across the wall. This should 
not be a bad approximation when there are many light species of particles 
in the system. One obtains for the f-dependent force on the wall, a result 
roughly suppressed by a factor of order \J~Cj~8- For a Higgs mass of 35 GeV 
we obtain a velocity ~ 0.2. 

In this thick wall case, the density enhancement described earlier could be 
extremely important. Here one expects the extra particles to be distributed 
more or less uniformly over the wall. This leads to an extra force AF ~ 
fnvr Ap /5 where Ap denotes the internal energy difference on the two sides 
of the wall. As an estimate of /, we can take the ratio of the equilibrium 
densities on the two sides of the wall; for top quarks, this gives a number 
of order 5%. This effect seems to be comparable to that of the preceding 
paragraph. 

7. In this paper we have outlined a program for treating the phase 
transition in weakly coupled theories with scalar fields. From the standpoint 
of electroweak baryogenesis, it is important to find models consistent with 
current experimental bounds in which the phase transition is strongly first 
order. It would be interesting to find models where the bubble wall is thin; 
in such theories || electroweak baryogenesis can be very efficient. However, 
our work suggests that in many models, the wall will be slow and thick, and 
adiabatic analyses of the type of Ref. will be relevant. In these cases, 
one may just barely be able to produce the observed asymmetry. Further im- 
provements in both the theory of the phase transition and that of electroweak 
baryogenesis will be necessary to completely settle these questions. 

We are grateful to L. Susskind, L. McLerran, R. Kallosh, N. Turok and 
M. Gleiser for discussions. 
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